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We study the static susceptibilities for charge and spin sectors in paramagnetic states for 
Ca2-iSr l; Ru04 in 0.5 < x < 2 within random phase approximation on the basis of an effective 
Ru t2 g orbital Hubbard model. We find that several modes of spin fluctuation around q = (0, 0) 
and q ~ (0.7977T, 0) are strongly enhanced for the model of x — 0.5. This enhancement arises from 
the increase of the corresponding susceptibilities for the d xy orbital due to the rotation-induced 
modifications of the electronic structure for this orbital (i.e., the flattening of the bandwidth and 
the increase of the density of state near the Fermi level). We also find that the ferromagnetic spin 
fluctuation becomes stronger for a special model than for the model of x = 0.5, while the compe- 
tition between the modes of spin fluctuation at q — (0, 0) and around q ~ (n, 0) is weaker for the 
special model; in this special model, the van Hove singularity (vHs) for the d xy orbital is located 
on the Fermi level. These results indicate that the location of the vHs for the d xy orbital, which 
is controlled by substitution of Ca for Sr, is a parameter to control this competition. We propose 
that the spin fluctuations for the d xy orbital around q = (0, 0) and q ~ (n, 0) play an important 
role in the electronic states around x — 0.5 other than the criticality approaching the usual Mott 
transition where all electrons are localized. 

PACS numbers: 71.27.+a, 74.70.Pq 



I. INTRODUCTION 

Ca 2 - a; Sr a ;Ru04 is one of the strongly correlated elec- 
tron systems with orbital degrees of freedom.EE This 
alloy has a quasi-two-dimensional crystalline structure, 
and the Ru t2 g orbitals and the lattice distortions in- 
duced by substitution of Ca for Sr play important roles 
in determining the electronic structures. Actually, this 
alloy shows a spin-triplet superconductivityJ^El a heavy 
fermion (HF) behavior,^ and a metal-insulator transi- 
tion,^ depending on the Sr concentration. 

To clarify the origin of the HF behavior around x = 
0.5, the present authors studied the electronic states for 
Ca2- 2: Sr 2: Ru04 in 0.5 < x < 2 within the Gutzwiller 
approximation on the basis of an effective Ru t2 g or- 
bital Hubbard modelP this effective model takes account 
of the change of the dp hybridizations due to the rota- 
tion of RuOg octahedra, which is induced in the range 
of 0.5 < x < 1.5. We proposed that this HF behavior 
can be qualitatively understood as the cooperative effect 
between moderately strong electron correlation and the 
orbital-dependent modification of the electronic struc- 
tures due to the rotation of RuC>6 octahedraP In that 
study however, we neglected the short-range correla- 
tions for charge and spin sectors, which will play a non- 
negligible role in discussing the electronic states. It is 
thus needed to study the effects of the short-range corre- 
lations on the electronic structures. 

In this paper, we study the static susceptibilities for 
charge and spin sectors in paramagnetic (PM) states for 
Ca2- :E Sr 2 ,Ru04 in 0.5 < x < 2 within random phase 
approximation (RPA) on the basis of the effective Ru t^g 
orbital Hubbard modelP In particular, we analyze these 
static susceptibilities for the models of x = 2 and 0.5, 



and a special model to analyze effects of the rotation of 
RuOg octahedra, the van Hove singularity (vHs) for the 
d xy orbital, and the Hund's rule coupling. In this special 
model, the vHs for the d xy orbital is located on the Fermi 
level. 

The paper is organized as follows. In Sec. II, we 
first show the effective Ru t2 g orbital Hubbard model 
for Ca 2 _ 2; Sr a .Ru04 in both the absence and the presence 
of the rotation of RuOg octahedra. Next, we introduce 
a generalized susceptibility for a multiorbital system and 
susceptibilities for charge and spin sectors. Then, we ex- 
plain the method to determine a primary instability for a 
multiorbital system. In Sec. Ill, we first show the static 
susceptibilities and the maximum eigenvalues of static 
susceptibility for the models of x = 2 and 0.5 to analyze 
the effects of the rotation of RuOg octahedra. We next 
compare the results for the model of x = 0.5 with those 
for the special model to analyze the effects of the vHs for 
the d xy orbital. Then, we show the dependence of the 
static susceptibilities on the Hund's rule coupling for the 
models of x = 2 and 0.5, and the special model. In Sec. 
IV, we first address the mass enhancement due to the en- 
hanced fluctuations for the models of x = 2 and 0.5. In 
addition, we compare our results with previous theoret- 
ical and experimental works about magnetic properties. 
The paper concludes with a summary of our results in 
Sec. V. 



II. FORMULATION 

In the following, we set h = /xb = &b — 1, label the 
d X z, d yz , and d xy orbitals as 1, 2, and 3, and choose the 
coordinates, x, y, and z, in the directions of the Ru-0 
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FIG. 1: FSs for the models of (a) x = 0.5 and (b) 2. The 
dashed line in Fig. (a) represents the folded BZ. The arrows 
in Figs, (a) and (b) correspond to the momentum lines used 
in the following figures. 



bonds at <j> = 0° (i.e., 1.5 < x < 2). Here, <fi is an angle 
of the rotation of RuOe octahedra. 

To study the electronic states for Ca2- x Sr a; Ru04 in 
0.5 < x < 2, we consider an effective Ru t 2g orbital 
Hubbard model which takes account of the change of the 
dp hybridizations between the Ru Ad and O 2p orbitals 
due to the rotation of RuOk octahedral 



H = Hn 



int. • 



(1) 



In the presence of the rotation of RuOg octahedra, the 
alternation of the direction of the rotation in the two- 
dimensional square lattice leads to a unit cell doubled, 
and the Brillouin zone (BZ) is folded [see Fig. [I] (a)]. 

The noninteracting Hamiltonian in the absence of the 
rotation of RuOg octahedra is given by 



#0=E E E t ab (k,o° 

k a,b=l s="f,4- 



^kas^kbs 



(2) 



where e a &(fc, 0°) denotes the energy dispersions for cj) — 
0°. These energy dispersions measuring from the chemi- 
cal potential, (j,, are given by 



+ cos k v ) — 4i4 cos k x cos k y — 



(3) 
(4) 
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en(fc, 0°) = —2t\ cos k x — 2t2 cos k y — fj,, 
ei 2 (fe,0°) = e 2 i(fe,0°) = -41/ sin k x sin k y , 
e 2 2(fc, 0°) = —2t 2 cos k x — 2ti cos k y — fi, 
£33(^,0°) = -2t 3 (cosk x 

e a b{k,0°) — otherwise 

In this study, we take account of the next-nearest- 
neighbor (n.n.n.) hopping integral between the d xz and 
d yz orbitals, which has been neglected in the previous 
study (Ref|H|. A* is determined so that the total occupa- 
tion number for the Ru ti g orbitals is equal to 4. 

For the model of x = 2, we set h = 0.675, t 2 = 0.09, 
<3 = 0.45, <4 = 0.18, and t' — 0.03 so as to reproduce the 
experimentally observed Fermi surfaces (FSs) P Here and 
throughout this paper, the energy unit is eV. Figure [l] 
(b) shows the FSs for this model. We see that the n.n.n. 
hopping integral between the d xz and d yz orbitals leads 
to the slight change of the topology of the FS around 



k ~ (0.677T, 0.677r) from that without this n.n.n. hopping 
integral (RefiS); there are no qualitative changes in the 
density of states (DOS), which is not shown here. 

In the presence of the rotation of the RuOg octahedra, 
the noninteracting Hamiltonian becomes 



= E E E 

k a,b=l l,l'=A,B s=f4 



E e aU k ^) C kals C kb- 



at 
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where the prime in the summation with respect to mo- 
mentum represents the restriction within the folded BZ 
for cj) 7^ 0°, and I and e l * b (k,(f>) denote the sublattice in- 
dex and the energy dispersions for 0^0°. These energy 
dispersions measuring from /i are given by 



e^{k,ct>) = e^{k,ct>) = \^ 



Mi 



efa (fc,0) = ejfi (k,4>) = -W sin k x sin h 



(9) 
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£3 3 A (fc, 4>) = - - A t2g - 4t 4 cos k x cos k y - fi, (11) 

e^ 8 (k,cf>) = — 2t\ cos 2 <fi cos k x — 2(t 2 — t\ sin 2 <j>) cos k 

(12) 

(13) 

(14) 
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ti sin 20(cos k x + cos k y ), 

-6f 2 B (fe,0), 



€-22 {k, 4 1 ) — ~ 2(ia — ti sin cj)) cosk x — 2t\ cos <pcosk y . 
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2^3 cos 3 2(/)(cos k x + cos k y 



+ 2t5 cos 2(j> sin 2 20(cos k x + cos k y ) 

— Atg cos 20 sin 2 20(cos k x — cos k y ), (16) 

e^(fc, 4>) = e* b B (k, 0) = otherwise, (17) 

£ f h s (fe,0) = e ^(fc,-0), (18) 

e3*(M) = 4 b B (k-ct>). (19) 



The detail of the derivation of these energy dispersions 
is described in Ref. HI In addition to the change of the 
hopping integrals, we have taken account of the effect 
of the rotation-induced hybridization of the d xy orbital 
to the d x 2_ y 2 orbital as the difference of the crystalline- 
electric-field (CEF) energies between the d xz /y Z and d X y 
orbitals, A t2g P 

For the model of x — 0.5, we set <p = 15°, At 2g = 0.39, 
and t5 = te = so as to reproduce the experimentally 
observed FSsP^Note that the value of A t2g in this study 
is different from that in the previous studies^^ due to 
the introduction of the n.n.n. hopping integral between 
the d xz and d yz orbitals. Figures [l] (a) and [2| s how the 
FSs and DOS for this model. We see from Fig. |2jthat the 
n.n.n. hopping integral between the d xz and d yz orbitals 
leads to a sharp decrease of the DOS for the d xz / yz orbital 
around an energy of —0.4 due to a small gap opening in 
the bands for these orbitals. This result is consistent 
with that obtained in the density-functional calculation 
for x = 0.5P 

In addition, to discuss the role of the vHs for the d xy 
orbital, we consider the special model, for which we set 
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FIG. 2: (Color online) DOS for the model of x = 0.5. The 
dashed black line represents the chemical potential. 



The noninteracting susceptibility is 
XT («)=~EE GT (k + *)< (0) (*). (25) 

Here, k = (fc, iuj m ) with cj,„ being a fermionic Matsubara 

frequency, and G l ^°\k) is the noninteracting Green's 
function: 



G l T(k)=UUk) al ; a - 



1 



lU) r , 



€ a (fc, 



-(Ul) a;W , (26) 



where (Uk)ai-.a is a unitary matrix to diagonalize Hq: 

Ckals — J2 a {Uk)al;aCkas and 

= E E (^)a;a/eL i ;(fc,0)(C/ fc ) 6 /'; a . (27) 
a,b=l Z,i'=A,B 



e Q (fc, 



^ = 15°, A t2s = 0.27, and = t e = 0; similarly to the 
model of x — 0.5, the value of A t2s is different from that 
in the previous studies (Refs. 151 and ITT]) . In this special 
model, the vHs for the d xy orbital is located on the Fermi 
level, as denoted in Sec. I. 

We assume that the interacting Hamiltonian is given 
by the standard on-site interactions: 



U 2~2J2 njaffijai + C/'EE n ja n jb 

j a j a>b 

-4EE ( 2 *ja ' Sjb + \fl ja fl jb ) 
j a>b 

+ J 'EE d ]at d ]al d 3blC jbt , 
j a>b 



where n,- 



E S , S ' 4as a s S ' d 3as' and s ja 



Es flj as 

|E s ,s' Cjas cr ss>Cjas> with a° s , and cr ss > being the unit 
matrix and the Pauli matrices. 

To analyze fluctuations for charge and spin sectors for a 
multiorbital system, we define a generalized susceptibility 
at temperature T as 

Xas\ ,bs2 ;CS3 ,ds4 

= E (^rc{ lalsi (T)c k > +qb i S2 (T)cl cl ,^C k ^ qdVsi ), (21) 
k,k' 

and the corresponding Fourier coefficient is 

rT- 1 ^ u , 

XaSl ,f)S2XS3 ,d,S4 (*?) / XaSi,bS2]CS3,dS4 (*7? ^~) - 

Jo 

(22) 

Here, T T is a chronological operator for imaginary time 
t, c k aia(T) = e rH c kala e- TH , and q = (g,iQ„) with Q„ 
being a bosonic Matsubara frequency. In a PM state, the 
generalized susceptibility satisfies 



(9), 
(?)■ 



(23) 
(24) 



Let us introduce susceptibilities for charge and spin 
sectors as follows: 



XabSKq) = x"t,bT;cT,dt(9) + Xai.b^ct.dfC?). ( 28 ) 

XabfJil) = Xal.bfjct.rftCe) " xfi,fc*;eMt(?)- ( 29 ) 

These susceptibilities characterize the fluctuations in a 
PM state since the spin-orbit interaction is neglected. 
For example, the fluctuation for spin degrees of freedom 
(i.e., spin fluctuation) is characterized by 



x s {q)=l E E xLfbii), 

L l,V=A,B a,b=l 



(30) 



(20) where 1/2 is the normalization constant for the summa- 



tion with respect to the sublattice indi ces. 

Xabfd\<l) .and x2'(!) in Ec l s ; B and Q are de- 
tcrmincd within the RPA by using the following equa- 
tionsPHUl 

l"l'(C) 



Y U ' {G) (a) - y ll ' {0) (a) 

Xabcd \1) — Xabcd W 



Y U ' iS) (a) - Y U ' {0) (a) 
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(31) 

rr y ""(») f(;lr s y l " l ' {s) (a) 

2—i 2-i Aaba'b'W 1 a'V c! d' A c < ' d< ' cd W> 
{a'} I" 

(32) 



where E{ a '} = E a ' b> c' d> an d the bare vertex interac- 
tions are 
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respectively. These bare vertex interactions are indepen- 
dent of the sublattice indices since the interacting Hamil- 
tonian (20) is diagonal in terms of these indices. 



Since a phase transition is characterized by a diver- 
gence of static susceptibility, we calculate the eigenvalues 
of the static susceptibility by solving 



E E Xa'bidHq, 0)u^l. v (q) = \°(q)u%}.M, (35) Q = t) approach q = (tt, tt) as x is changed from 2 



the cases of x — 2 and 0.5.] We find that the noninter- 
acting susceptibility for the d xy orbital around q = (0, 0) 
is larger for the model of x — 0.5 than for the model of 
x — 2. This difference is due to the increase of the DOS 
for the d xy orbital near the Fermi level in the presence of 
the rotation of RuC>6 octahedraP We also find that the 
incommensurate (IC) peaks for the d xz / yz orbital around 



■d I' 



EEx2(9,0» ;1/ (g) = Xl(q)u^. iV (q). (36) 

c.d V 

Here, X^ S (q) and w^.^(q) are the eigenvalue of the 
static susceptibility for a charge/spin sector and the cor- 
responding eigenvector. We analyze the instability from 
the maximum eigenvalues!^ 
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Finally, to address the role of each Ru ti g orbital in the 
fluctuations, we introduce the susceptibilities averaged 
with respect to the sublattice indices: 
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^abct 
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r, Xabcd W- 
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(39) 
(40) 
(41) 



III. RESULTS 

In the following calculations, we set T = 0.02, J' = Jh, 
and V = U — 2Jn, and use the values of U and Jh as 
parameters. To calculate the static susceptibilities, we di- 
vide the BZ into 128 x 128 meshes and take 1024 fermionic 
Matsubara frequencies for using the fast Fourier transfor- 
mation. 



A. Effects of the rotation of RuOg octahedra 

To analyze the effects of the rotation of RuOg octa- 
hedra on the static susceptibilities, we compare the re- 
sults for the model of x — 2 (i.e., <j> = 0°) and 0.5 (i.e., 
0=15°). 

First, the momentum dependences of the noninteract- 
ing susceptibility for these models are shown in Figs. [3] 
(a) and (b). [Note that the momentum lines in both fig- 
ures are along the arrows shown in Fig. [I] (b) to compare 



to 0.5; simultaneously, the IC peak for the d xy orbital 
around q = (tt,0) approach q = (tt,0). These shifts are 
understood as follows. The rotation-induced hybridiza- 
tion of the d xy orbital to the d x i_ y i orbital leads to the 
downward shift of the energy level for the d xy orbital. 8 
Correspondingly, the energy level for the d xz / yz orbital 
goes up compared with that for the d xy orbital. The 



upward shift of the band for the d x 



orbital around 



q = (tt,tt) crossing the Fermi level leads to an increase 
of the nesting vector towards q = (it, it). 

Next, we show the momentum dependences of the max- 
imum eigenvalue of static susceptibility in the RPA at 
Jh = U/6 for the models of x — 2 and 0.5. Figure [4] (a) 
shows the momentum dependences of A^ lax (<7) _1 for the 
model of x = 2. We find that A^ lax (q)~ 1 at q ~ (tt, 0.5tt) 
touches zero at U = U c = 0.975, and that there are two 
main peaks at q ~ (tt, 0.5it) and (0.2357T, 0.2357r) and two 
secondary peaks at q ~ (tt, 0.6567t) and (0.6887T, O.68871-). 
Note that the enhancement of these main peaks is in- 
consistent with the experimental result of an inelas- 
tic neutron measurement,^ in which the enhancement 
of the IC antiferromagnetic (AF) spin fluctuation with 
q ~ (0.67T, O.671-) is observed. However, as we will show 
below, our model can explain this experimental result. 
We also find that the increase of U leads to the sup- 
pression of the values of A^ ax (q) _1 (not shown here). It 
should be noted that the small value of U c is due to the 
overestimation of the fluctuations in the RPA. This small 
value of U c will be unrealistic since the value of U c for 
transition metal oxides will be of the order of 2 — 3 eV. 
The similar problem remains for the model of x = 0.5 
and the special model, as we will show below. In Sec. 
IV, we will address this problem and effects of the terms 
neglected in the RPA. 

The momentum dependences of A^ lax (q) _1 for the 
model of x = 0.5 are shown in Fig. [4] (b); the inset 
represents A^ ax (q) _1 at U = U c = 0.751 in an expanded 
scale. [Hereafter, the momentum lines for 0^0° are 
along the arrows shown in Fig. [I] (a).] We find that 
Amax(9) -1 at g = (0,0) touches zero at U = U c = 0.751, 
which is smaller than for the model of a; = 2, and that 
there are two main peaks at q = (0, 0) and q ~ (0.7977r) 
and a secondary peak at q ~ (tt, 0.1257r); all these peaks 
are less sharp than the peaks for the model of x = 2. 
As we will discuss in Sec. IV B, the enhancement of the 
main peak at q = (0, 0) is consistent with several exper- 
imentsPHini \y e also find that the values of A^ ax (q) _1 
(not shown here) suppress as U increases. These results 
are qualitatively same to those obtained in the previous 
studyP where the n.n.n. hopping integral between the 
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FIG. 4: (Color online) Momentum dependences of A^ax 

special model. The insets in (b) and (c) show the momentum dependences of A 
the model of x = 0.5 and the special model, respectively. 



(q) 1 in the RPA for the models of (a) x = 2 and (b) 0.5, and (c) the 

(q) _1 at U — U c in an expanded scale for 



s 

max V 



d xz and d yz orbitals has been neglected. 

Moreover, Figs. [5~| (a) and (b) show the momentum 
( sv — 

dependences of \aaaa (q, 0) in the RPA at U = 0.8U C for 
the models of x — 2 and 0.5. For the model of x = 2, 
we see from Figs. [4] (a) and [5] (a) that the two main 
peaks in A^ ax (q) _1 arise from the corresponding fluctu- 
ations for the d xy orbital, and that the secondary peaks in 
Af nax (q)" 1 at q ~ (7r,0.6567r) and (0.688tt, 0.688tt) arise 
from the corresponding fluctuation for the d yz orbital and 
the combined fluctuation of the d xz / yz and d xy orbitals. 
These results suggest that the fluctuations for a spin sec- 
tor not only for the d xy orbital but also for d xz / yz orbital 
play an important role for x = 2. 

For the model of x — 0.5, we see from Figs.[|(b),|5](b), 
andpka) that the two main peaks and the secondary peak 
in A^ ax (q) _1 arise from the corresponding fluctuations 
for the d xy orbital, and that the contribution from the 
d xz / yz orbital is much smaller than for the d xy orbital. 

Note that the value of Xaaaa\q, 0) for the d xy orbital is 
about 10 times larger than for the d xz / yz orbital. These 
results suggest that the rotation of RuC>6 octahedra leads 
to the enhancement of several modes of spin fluctuation 
for the d xy orbital around q = (0,0) and q ~ (7T,0), 
and that the fluctuation for a spin sector plays a more 
important role for x = 0.5 than for the d xz / yz orbital. 
Finally, Figs. [7] (a) and (b) show the momentum de- 



pendences of x S (?:0) m the RPA at Jh = U/6 for U — 
0.8U C and 0AU c for the models of x = 2 and 0.5. For the 
model of x = 2, we find from Fig. [7] (a) that x S (<Zi 0) at 
q ~ (vr, 0.5tt), (tt, 0.656tt), and (0.688tt, 0.688tt) are domi- 
nant at U = 0.8U C , while x s (<7, 0) at q ~ (0.688tt, 0.688tt) 
becomes largest at U = 0AU c ; the main peak in 
A^ ax (q) -1 at q ~ (0.235-7T, 0.2357r) is less important at 
U = 0.8U C and 0.4J7 C . These results suggest that the 
experimentally observed enhancemenlP^ of the 1C AF 
spin fluctuation at q ~ (0.6-7T, 0.67r) can be understood 
if Sr2RuC>4 is not located on the vicinity of the magnetic 
order. As we will show in Sec. HI C, the Hund's rule 
coupling also plays the important role in enhancing this 
IC AF spin fluctuation. 

For the model of x = 0.5, we find from Fig. [7] (b) that 
X S (<7, 0) at q = (0, 0) and q ~ (tt, 0.1257r) are dominant at 
U = 0.8U C , while x S (<?,0) at q = (0,0), q - (0.797tt,0), 
and (71", 0.1257r) are nearly same at U = 0.4LV,,. Namely, 
almost all contributions arise from the corresponding 
fluctuations for the d xy orbital. We also find that the 
values of x S (<Zj0) along (0,0) — > (n,0) are larger than 
along (0,0) — > (7r/2,7r/2). This result indicates that the 
spin fluctuations along (0, 0) — > (n, 0) are strongly en- 
hanced for x = 0.5. This result i s con sistent with several 
inelastic neutron measurements GH2H] as we will address 
in Sec. IV B. These results suggest that the spin fluctua- 
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(a)x = 2 (b)x = 0.5 (c)vHsonE F 




FIG. 5: (Color online) Momentum dependences of Xaaaa(q, 0) in the RPA for the models of (a) x — 2 and (b) 0.5, and (c) the 
special model. 
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FIG. 6: (Color online) Momentum dependences of 
Xwaa (q, 0) for the d xy orbital in the RPA for (a) the model 
of 0.5 and (b) the special model. 

tions for the d xy orbital along (0, 0) — > (tt, 0) play a very 
important role for x = 0.5. 



B. Effects of the vHs for the d xy orbital 

To analyze the effects of the vHs for the d xy orbital on 
the static susceptibilities, we compare the results for the 
model of x = 0.5 with those for the special model, where 
the vHs is located on the Fermi level; the difference in 
these models is only the value of A t2 . 

First, the momentum dependence of the noninteracting 
susceptibility for the special model is shown in Fig. [3](c). 
Comparing this figure with Fig. [3] (b), we see that the 
noninteracting susceptibility for the d xy orbital around 
q = (0, 0) is larger for the special model than for the 
model of x — 0.5. This difference is due to the increase 
of the DOS for the d xy orbital near the Fermi level in 
the special model. We also see that the shifts of the IC 
peak for the d xz / yz and d xy orbitals, which are observed 
in the model of x = 0.5, are smaller in the special model. 
This difference arises from the smaller value of A*, for 
the special model. 

Next, Fig. [4] (c) shows the momentum dependences 
of Af nax (q)^ 1 in the RPA at J H = U/6 for the special 
model; the inset represents A^ ax (q) _1 at U = U c = 0.701 
in an expanded scale. Comparing Fig. [I] (c) with Fig. [4] 
(b), we find that the value of U c , where A^ ax (q) _1 at 
q = (0, 0) touches zero, is smaller for the special model 
than for the model of x = 0.5, and that there are a 



main peak at q — (0,0) and a secondary peak at q ~ 
(0.7037T, 0); the latter peak corresponds to the peak at 
q ~ (0.797-7T, 0) for the model of x = 0.5. In addition, 
we see that all the peaks in A|j ax (q) _1 for the special 
model are slightly sharper than for the model of x — 
0.5, and that the competition between the modes of spin 
fluctuation around q — (0, 0) and q ~ (tt, 0) is weaker for 
the special model. The latter is due to the larger increase 
of the noninteracting susceptibilities for the d xy orbital 
around q = (0, 0) compared with those around q ~ (tt, 0) 
for the special model than for the model of a; = 0.5. 
These results indicate that the location of the vHs for 
the d xy orbital, which is controlled by substitution of 
Ca for Sr, is the parameter to control this competition. 
Note that the momentum dependences of A£ ax (q) _1 and 
Xa aaa (q,0) (not shown) are qualitatively same to those 
obtained for the model of x — 0.5. 

Finally, Figs. [5] (c), [6] (b), and [7] (c) show the mo- 
mentum dependences of Xaaaa(q,0), X3333 S) ('?> )> and 
X S (<7,0) in the RPA for the special model. Comparing 
these figures with Figs. [5] (b), [6] (a), and [7] (b), we find 
the qualitatively same results to those for the model of 
x = 0.5. Namely, the enhanced modes in A^ lax (q) _1 arise 
from the corresponding fluctuations for the d xy orbital, 

and the dominant contributions to Xauaa(<l,ti) arise from 
those for the d xy orbital around q = (0, 0) and q ~ (tt, 0). 
These results suggest that the spin fluctuation for the d xy 
orbital is more important in the electronic states around 
x = 0.5 than for the d xz i yz orbital. 

C. Effects of the Hund's rule coupling 

To analyze the effects of the Hund's rule coupling on 
the static susceptibilities, we compare the results in the 
RPA at J H = U/6 with those at J H = and U/A. 

First, we focus on the Jh dependence of the static sus- 
ceptibilities for the model of a; = 2. We find that the 
dominant wave vectors of these susceptibilities at Jh = 
and U/A (not shown) are same to those at Jh = U/6, 
and that the increase of Jh/U leads to the decrease of 
U c . The latter is a typical behavior for a multiorbital 
system since the increase of Ju/U leads to an enhance- 
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FIG. 8: (Color online) Momentum dependences of Xaaaa(q, 0) 
in the RPA for the models of (a) x — 2 and (b) 0.5. 



merit of spin fluctuations. We also find from Fig. [8] 

(a) that the static susceptibility for a charge sector is 
enhanced only at Jjj = as U increases. The sim- 
ilar enhancement has been reported in a RPA analy- 
sis for x = 2pH where the n.n.n. hopping integral be- 
tween the d xz and d yz orbitals has been neglected. Fur- 
thermore, we see from Fig. [91(a) that x s (<7,0) around 
q ~ (0.688-7T, O.68871-) and (7r,0.6567r) are strongly en- 
hanced as J-r/U increases; these modes correspond to the 
secondary peaks in A^ ax (q) _1 at Jh = U/6. The nature 
of these enhancements can be understood as follows!^ 
the Hund's rule coupling, which is one of the interorbital 
Coulomb interactions, more affects the orbitals having a 
finite hybridization. 

Next, we turn to the Jh dependence of the static sus- 
ceptibilities for the model of x = 0.5. Similarly to the 
case of x = 2, we find that the dominant wave vectors of 
the static susceptibilities at Jh = and U/4 (not shown) 
are same to those at Jh = U/6, that the increase of 
Je_/U leads to the decrease of U c , and that the increase 
of U leads to the enhancement of the static susceptibil- 
ity for a charge sector only at Jh = [see Fig. [8] (b)]. 
In contrast to the case of a; = 2, we find from Fig. [9] 

(b) that the increase of Jh/U leads to a nearly uniform 
enhancement of x S (<7;0); more precisely, the enhance- 
ment of x s (q, 0) around q ~ (it, 0) is slightly larger than 
around q = (0,0). This difference is due to the much 
larger contribution from the d xy orbital than from the 
d-xz/yz orbital for the model of x = 0.5. These results in- 



dicate that the static susceptibilities for a spin sector for 
the d xy orbital are very important for x = 0.5 regardless 
of the value of Jr/U. 

Before showing the results for the special model, let us 
remark on the effect of the rotation-induced hybridiza- 
tion of the d xy orbital to the d x 2_ y 2 orbital on the Jh 
dependence of x s (<7,0) for the model of x = 0.5. In this 
study, the effect of this hybridization is approximately 
taken account as the difference of the CEF energies be- 
tween the d xz / yz and d xy orbitals [see Eqs. (|9| and ( 11 )]. 
In principle, if this hybridization is fully taken into ac- 
count, the Jh dependence of x s (g, 0) will be modified. 
However, according to the density-functional calculation 
for x = 0.5p^ the d x 2_ y 2 orbital has little weights to 
the bands near the Fermi level for x = 0.5. This re- 
sult indicates that the values of the static susceptibility 
for the d x 2_ y 2 are much small compared with those for 
the tig orbitals. Thus, the Jh dependence of X S (q, 0) 
does not qualitatively change from the present result if 
this rotation-induced hybridization is fully taken into ac- 
count. 

Finally, we focus on the Jh dependence of the static 
susceptibilities for the special model. We find the similar 
Jh dependences of the value of U c and the static suscep- 
tibilities to those for the model of x = 0.5 [e.g., see Fig. 
[91(c)]. Therefore, we conclude that the fluctuations for a 
spin sector for the d xy orbital play a very important role 
in the electronic states around x = 0.5. 



IV. DISCUSSION 

We first address the mass enhancement for the models 
of x — 2 and 0.5. In this study, we find that the value 
of U c is smaller for the model of x = 0.5 than for the 
model of x = 2, and that the peaks in A^ ax (q) _1 are 
less sharp for the model of x = 0.5 [see Figs. [4] (a) and 
(b)]. In general, the effective mass at a certain value of 
U becomes large if the value of U c is smallP^ID j n a j_ 
dition, the effective mass becomes larger for a less sharp 
peak in the susceptibility, even if the value of the peak 
is samepSHUl this arises from the product of the suscep- 
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FIG. 9: (Color online) Momentum dependences of X S (<7iO) m the RPA for the models of (a) x = 2 and (b) x = 0.5, and (c) 
the special model. 



tibility and the summation with respect to q in the free 
energy. Combining these with our results, we conclude 
that the mass enhancement for the model of x — 0.5 will 
be much larger than for the model of x — 2. Namely, we 
think that the present results can qualitatively explain 
the experimentally observed mass enhancement towards 
x = 0.5 (RefjS]); the actual calculation about the mass 
enhancement is a remaining future work. 



A. Comparison with previous theoretical studies 

We first compare our results with the previous 
theoretical studies about the magnetic properties for 
Ca2- :! ;Sr 2: Ru04 with x = 2. The mean-field calcula- 
tion^! for x = 2 has found that the dominant contri- 
bution to the static susceptibility arises from that for 
the d xy orbital. This calculation has also shown that 
the increase of Jn/U leads to the enhancement of the 
IC AF spin fluctuation around q ~ (0.677T, 0.677r), and 
that this IC AF spin fluctuation is strongest only for 
Jh > U/2 at U = 1.25t 3 ~ 0.56. Note that this 
mode of the enhanced IC AF spin fluctuation is differ- 
ent from q ~ (0.6887T, 0.6887r), which is enhanced in our 
study, due to the difference in the values of the hop- 
ping integrals. These results suggest that a moderately 
large Hund's rule coupling plays an important role for 
x = 2. This is consistent with the previous dynam- 
ical mean-field theory^ for x = 2. Similarly to this 
mean- field studyp^ we find that the increase of Jb_/U 
leads to the enhancement of the IC AF spin fluctuations 
around q ~ (0.6887T, 0.6887r) and (tt, 0.6567r); the for- 
mer IC AF spin fluctuation is strongest for Jh > U/6 
at U = 0.8U C ~ 0.78. Since the enhancement of the 
IC AF spin fluctuation at q ~ (0.67T, 0.67T, 0) has been 
experimentally observed,^ we think that a set of the pa- 
rameters used in the present study is more realistic than 
in this previous mean- field study (Rcf.22). 

Next, we address the physical meaning of the re- 
sults obtained in the RPA, and compare our results 
with the studies in more elaborated treatments. The 
RPA is a mean-field type approximation and neglects 
the mode-mode coupling for fluctuations. In general, 



the mode-mode coupling plays an important role in dis- 
cussing the electronic states near a quantum critical point 
(QCP)P^] it weakens the instability obtained in the 
RPA. Actually, the previous study^for x = 2 in fluctua- 
tion exchange (FLEX) approximation has observed much 
larger U c than in the RPA; the FLEX approximation is 
partially taken account of the mode-mode coupling. We 
thus expect that the values of U c will be more realistic 
than those obtained in this paper. 

In addition to this, the mode-mode coupling can 
change the primary instability from that obtained in the 
RPA. Thus, the results about the primary instability may 
change when the mode-mode coupling is taken into ac- 
count. Actually, the previous study^ for x = 2 in the 
FLEX approximation found that the IC AF spin fluctua- 
tion with q ~ (0.677T, 0.677r) is very strong, while the spin 
fluctuations with q ~ (tt, 0.5tt) and (0.2357T, 0.2357r) are 
strongly enhanced in the RPA; the reason why this mode 
of the enhanced IC AF spin fluctuation in the FLEX is 
different from that in our study is same to that for a 
case of the difference between the previous mean-field 
study^l and ours. However, we believe that the com- 
petition between the modes of spin fluctuation around 
q = (0, 0) and q ~ (tt, 0) for the model of x = 0.5, 
which is obtained in the present study, will not change 
even in a more elaborated treatment since the momen- 
tum dependence of A^ lax (q) _1 is very flat; in such a case, 
it will be difficult that the mode-mode coupling favors 
a specific mode. Actually, our preliminary calculation^ 
in the FLEX approximation observes the similar com- 
petition of spin fluctuations around these modes for the 
model of x = 0.5. 

There are several density-functional calculations to 
study the ground states for 0.5 < x < 2. Among them, 
a density-functional calculation within local-density ap- 
proximation (LDA) has found that the rotation of RuOe 
octahedra enhances the ferromagnetic (FM) instability.^ 
Another density-functional calculation^ for x = 0.5 
within the LDA has predicted the nesting instability with 
q ~ (0.297T, 0.297T, 0). The authors have proposed that 
this nesting instability will lead to the mass enhance- 
ment around x = 0.5; this nesting vector is related to a 
new FS induced by the rotation-induced hybridization of 
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the d xy orbital to the d x 2. y 2 orbital. In contrast, another 
density-functional calculation for x — 0.5 within local 
spin-density approximation has found a disappearance 
of this nesting vector in the presence of the spin-orbit 
interaction, which is equal to 0.167 eV.^ In our study, 
this new FS does not appear since we have not treated 
directly the rotation-induced hybridization of the d xy or- 
bital to the d x 2_ y 2 orbital; instead, we take account of the 
effect of this hybridization as the difference of the CEF 
energies between the d xz / yz and d xy orbitals, as denoted 
in Sec. II. It is thus necessary to study the effects of 
this new FS on the electronic structures around x = 0.5. 
However, we think that the presence of this new FS will 
not lead to the drastic changes from the present results 
since this rotation-induced FS is small and the DOS for 
each Ru t2 g orbital obtained in the present study is al- 
most same to that obtained in these density-functional 
studies (Refs. [[2 and [35]). 

Finally, we compare the obtained results with the pre- 
vious studj0 in the Gutzwiller approximation. This pre- 
vious study calculated the total renormalization factor 
of the kinetic energy for the Ru t2 g orbitals, which is in- 
versely proportional to the mass enhancement, and found 
that the inverse of the total renormalization factor is 
largest for the model of x = 0.5. Namely, these results 
can reproduce the experimentally observed tendency of 
the effective mass in 0.5 < x < 2 (Ref. [SJ> . These re- 
sults suggest that the criticality approaching the usual 
Mott transition plays the important role in enhancing 
the effective mass towards x = 0.5; in this usual Mott 
insulator, the occupation numbers for the d xz / yz and d xy 
orbitals are 1 and 2, respectively. On the other hand, in 
the present study, we find that several modes of spin 
fluctuation for the d xy orbital around q = (0, 0) and 
q ~ (0.797-7T, 0) are strongly enhanced for the model of 
x = 0.5. As discussed above, this enhancement of sev- 
eral modes will lead to the larger mass enhancement for 
x = 0.5 than for x = 2. This mechanism of the mass 
enhancement can coexist with the proposal in the pre- 
vious study^ in the Gutzwiller approximation since the 
mass enhancement due to the criticality approaching the 
usual Mott transition arises mainly from that for the d xz 
and d yz orbitals. In other words, not only the critical- 
ity approaching the usual Mott transition but also the 
spin fluctuation for the d xy orbital around q = (0, 0) and 
q ~ (7r, 0) play an important role in enhancing the effec- 
tive mass towards x — 0.5. 



B. Comparison with experimental results 

We first address the correspondence of our results with 
several experimental results at x — 2. In addition to the 
inelastic neutron measurement,^ our results are consis- 
tent with the nuclear magnetic resonance (NMR) mea- 
surement This experiment has found that not only the 
IC AF spin fluctuation at q ~ (0.67T, 0.67r) but also the 
FM spin fluctuation plays the non-negligible role in en- 



hancing the nuclear spin-lattice relaxation rate at a Ru 
site. This result indicates that the FM spin fluctuation 
also plays an important role in discussing the magnetic 
properties for x — 2. Actually, we see that the spin fluc- 
tuation around q = (0, 0) is non-negligible, though the 
dominant components arise from the IC spin fluctuations 
around q ~ (ir,0.5ir), (0.235-7T, 0.2357r), (tt, 0.6567r), and 
(0.6887r,0.6887r). 

Next, we compare our results with several experimen- 
tal results in the presence of the rotation of RuO@ oc- 
tahedra. A recent muon spin relaxation (/xSR) mea- 
surement!^ for Ca2_ K Sr2;Ru04 with x — 0.5 and 1.5 
and Sr2Rui_j,Ti a 04 with y = 0.09 has shown that 
the muon relaxation rates rapidly increase as tempera- 
ture decreases, and that the zero-field ^SR time spectra 
are similar to those for dilute-alloy spin glasses. The 
authors have proposed the static spin glass order not 
only in Sr2Rui_ a Tij,04 but also in Ca2- K Sr a; Ru04 for 
0.2 < x < 1.6. However, all the previous experimental 
studies fo r Ca 2 _ x Sr 2; Ru0 4 in 0.5 < x < 2 (Refs. 1112161151 - 
I20)38|39]) have suggested that the ground states are the 
PM metals. In addition, the elastic neutron measurement 
for x = 1.5 by the authors for this /iSR measurement has 
observed the absence of the ordered moment for x = 1.5, 
while the /iSR measurement has estimated the ordered 
moment at 0.25/ie per Ru. Therefore, we think that the 
conclusion for this //SR measurement is doubtful, and 
that it is reasonable to analyze the static susceptibilities 
in the PM states to discuss the electronic states around 
x = 0.5. 

We now discuss the correspondence with several neu- 
tron measurements. The polarized neutron diffraction 
measurement!^ for x = 0.5 has observed the anisotropic 
spin density distribution at a Ru site, which is flatted 
along c axis. This result suggests that the main part 
of the spin density arises from that for the d xy orbital. 
Namely, the spin fluctuation for the d xy orbital is pre- 
dominant at x = 0.5. Actually, we find that the domi- 
nant contributions to the spin fluctuation arise from that 
for the d xy orbital, which are consistent with this exper- 
imental result. 

In addition to this neu tron measurement, the inelas- 
tic neutron measurements^^ f or x = q 52 have found 
that the FM scattering is maximum at an energy transfer 
of 0.4 meV, while the IC scatterings with q = (0.12, 0, 0) 
and (0.27, 0, 0) evolve at an energy transfer of ~ 2.5 meV. 
These results indicate that the FM spin fluctuation is 
predominant at lower energy transfer, while the IC spin 
fluctuations evolve as an ener gy tr ansfer increases. Fur- 
thermore, these measurements^^! have observed the in- 
tensity spreading in momentum space around the maxi- 
mum peak and the large amplitude of the real part of the 
dynamic susceptibility at the maximum peak. These re- 
sults indicate that several modes of spin fluctuation are 
more strongly enhanced for x — 0.62 than for x = 2, 
and that Ca2- 2: Sr 2 ,Ru04 with x — 0.62 is located nearer 
to the boundary of the magnetic instability. Therefore, 
all the results for the model of x = 0.5 and the special 
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model are qualitatively consistent with these experimen- 
tal results except the evolution of the IC scatterings; it 
is necessary to study the origin of the evolution of these 
IC scatterings. 

Our result for the model of x = 0.5 is also consistent 
with the NMR measuerementp^ in which the evolution 
of the FM spin fluctuation has been observed towards 
x = 0.5. This result of the NMR measurement is consis- 
tent with the dependence of the Wilson ration on the Sr 
concentration.^ 

We now briefly remark on the metamagnetic transition 
ob serve d for Ca2-xSr ;c Ru04 in the range of 0.2 < x < 
q 5)2141] j n present study, we focus on the electronic 
states for 0.5 < x < 2 only in the absence of the external 
magnetic field. It is necessary to study the effects of the 
external magnetic field on the short-range correlation; 
this issue is a remaining future work. 

Finally, we remark on the Curie- Weiss (CW) behav- 
ior in the spin susceptibility with a Curie constant corre- 
sponds to nearly 5 = 1/2 for 0.5 < x < 2PElThe origin of 
this CW behavior has not been clarified yet. In general, 
there are two mechanisms in emerging the CW-type tem- 
perature dependence in the spin susceptibility! 30 1 31 ! One 
is due to the formation of the localized moment, and the 
other is due to the drastic renormalization of the temper- 
ature dependence for the free energy by dynamic fluctua- 
tions. In the latter mechanism, the mode -mode coupling 
leads to this drastic renormalizationPHHUxhe emergence 
of the CW behavior in Ca2- a; Sr. E Ru04 will arise from 
the latter mechanism since all electro ns fo r the Ru tig 
orbitals are itinerant in 0.5 < x < 2 ) 10 l 42 l In addition, 
we find for the model of x = 0.5 that several modes of 
spin fluctuation for the d xy orbital around q = (0, 0) and 
q ~ (0.797-7T, 0) are strongly enhanced, and that the dom- 
inant contributions to the susceptibilities for a spin sector 
arise from those for the d xy orbital. Therefore, we think 
that the mode-mode coupling for the d xy orbital will play 
a more important role in discussing the spin susceptibil- 
ity around x = 0.5 than for the d xz / yz orbital. It is 
necessary to study the electronic states around x = 0.5 
on the basis of the microscopic theory taking account of 
the mode-mode coupling. 



V. SUMMARY 

To clarify the effects of the short-range correlations 
on the electronic structures for Ca2- a; Sr.j : Ru04 around 
x = 0.5, we studied the static susceptibilities for charge 
and spin sectors in the PM states within the RPA for the 
models of x = 2 and 0.5, and a special model. In par- 
ticular, we analyzed the effects of the rotation of RuOe 
octahedra, the vHs for the d xy orbital, and the Hund's 
rule coupling on these static susceptibilities. 

First, we analyzed the effects of the rotation of RuOe 
octahedra by comparing the results at Jh = U /6 for 
the models of x — 2 and 0.5. In the absence of in- 
teractions, we found that the rotation of RuOg octa- 



hedra leads to the increase of the noninteracting sus- 
ceptibility for the d xy orbital around q = (0, 0) and 
the shifts of the IC peak for the d xz / yz and d xy or- 
bitals towards q = (tt,tt) and (tt,0). For the model 
of x = 2 in the RPA, we found that A 1 ^ lax (q) _1 at 
q ~ (tt, 0.5tt) touches zero at U = U c = 0.975, and 
that there are two main peaks in A^ lax (<7) _1 at q ~ 
(tt, 0.57r) and (0.235-7T, 0.2357r) and two secondary peaks 
at q ~ (7r,0.6567r) and (0.6887T, 0.6887r); the two main 
peaks arise from the corresponding fluctuations for the 
d xy orbital, and the secondary peaks at q ~ (tt, 0.6567r) 
and (0.6887T, O.68871-) arise from the corresponding fluctu- 
ation for the d yz orbital and the combined fluctuation of 
the d xz / yz and d xy orbitals. We also found that x s (q, 0) 
with q ~ (7r,0.57r), (7r,0.6567r), and (0.688tt, 0.688tt) 
are dominant at U = 0-8U c , while x (<Zj0) with q ~ 
(0.688-7T, 0.688-7T) becomes largest at U = 0.4{/ c ; the main 
peak in A^ lax (q) _1 at q ~ (0.2357T, 0.2357r) is less impor- 
tant at U = 0.8U C and 0AU c . For the model of x = 0.5 in 
the RPA, we found that A^ lax (<7) _1 at q = (0, 0) touches 
zero at U = U c = 0.751, which is smaller than for the 
model of x = 2, and that there are two main peaks in 
Amax(9) 1 at q = (0,0) and q ~ (0.7977T, 0) and a sec- 
ondary peak at q ~ (it, 0.1257r); in contrast to the case 
of x = 2, all the peaks in A^ ax (<7) -1 arise from the cor- 
responding susceptibilities for the d xy orbital. We also 
found that x S (q> 0) at q = (0, 0) and q ~ (tt, 0.1257r) are 
dominant at U = 0.8U C , while x S (Q,Q) at q = (0,0), 
q ~ (0.7977T, 0), and (tt, 0.1257r) are nearly same at 
U = 0AU c ; the values of x s (q,0) along (0,0) (tt, 0) 
are larger than along (0,0) — > (tt/2, it/2). Furthermore, 
we found that the peaks in A^ ax (q) _1 for the model of 
x = 0.5 are less sharp than for the model of x — 2. These 
less sharp peaks and the smaller value of U c will lead to 
a larger effective mass for the model of x = 0.5 than for 
the model of x = 2. 

Second, to analyze the effects of the vHs for the d xy 
orbital, we compared the results at Jh = U/6 for the 
special model with those for the model of x = 0.5. In 
the absence of interactions, we found the larger increase 
of the noninteracting susceptibility for the d xy orbital 
and the smaller shifts of the IC peak for the d xz i yz and 
d xy orbitals towards q — (tt, tt) and (tt, 0) in the spe- 
cial model. The former is due to the larger increase of 
the DOS near the Fermi level for the d xy orbital, and 
the latter is due to the smaller value of At 2 . In the 
presence of interactions, we found that the value of U c , 
where A,^ lax (q)^ 1 at q = (0,0) touches zero, is smaller 
for the special model than for the model of x — 0.5, and 
that there are a main peak at q = (0, 0) and a secondary 
peak at q ~ (0.7037T, 0); the latter peak corresponds to 
the peak at q ~ (0.7977r,0) for the model of x = 0.5. 
We also found that all the peaks in A^ lax (q) _1 for the 
special model are slightly sharper than for the model of 
x = 0.5, and that the competition between the modes 
of spin fluctuation around q = (0, 0) and q ~ (tt, 0) is 
weaker for the special model. Furthermore, we found 
that all the enhanced modes in A^ ax (q) _1 arise from the 
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corresponding fluctuations for the d xy orbital, and that 
the dominant contributions to the static susceptibility for 
a spin sector arise from those for the d xy orbital around 
q = (0,0) and q ~ (n, 0); these results are qualitatively 
same to those for the model of x = 0.5. 

Third, we analyzed the Jh dependences of the static 
susceptibilities for the models of x — 2 and 0.5, and the 
special model. For all these models, we found that the 
dominant wave vectors of the static susceptibilities at 
Jh = and U/4 are same to those at Jh = U/6, that 
the increase of Jn/U leads to the decrease of U c , and 
that the static susceptibility for a charge sector is en- 
hanced only at Jh = 0. In addition to these, we found 
for the model of x = 2 that the increase of Jn/U enhances 
the IC AF spin fluctuations at q ~ (0.688tt, 0.688tt) and 
(ir, 0.6567r), and that the former IC AF spin fluctuation 
is strongest for J H > U/6 at U = 0.8U C ~ 0.78; these 
modes correspond to the secondary peaks in A,„ ax (<7) _1 
at Jh = U/6. The enhancement of the IC AF spin fluc- 
tuation at q ~ (0.6887T, O.68871") arises not only from the 
dxz/yz orbital but also from the d xy orbital. Although 
the similar results have been obtained in the previous 
mean- field studjE2l for x = 2, we think that a set of the 
parameters used in this study is more realistic. In con- 
trast to the case of x = 2, the increase of Jh/U leads to a 
nearly uniform enhancement of x s (q, 0) for the model of 
x = 0.5 and the special model. This difference is due to 
the much larger contribution from the d xy orbital than 



from the d xz / yz orbital in the presence of the rotation of 
RuOg octahedra. 

In summary, the analyses in the RPA reveal that the 
rotation of RuOg octahedra leads to the enhancement 
of several modes of spin fluctuation for the d xy orbital 
around q = (0,0) and q ~ (tt,0). This enhancement 
arises from the increase of the corresponding suscepti- 
bilities for the d xy orbital due to the rotation-induced 
modifications of the electronic structure for this orbital 
(i.e., the flattening of the bandwidth and the increase of 
the DOS near the Fermi level) . These analyses also reveal 
that the location of the vHs for the d xy orbital, which is 
controlled by substitution of Ca for Sr, is a parameter to 
control the competition between the modes of spin fluctu- 
ation for the d xy orbital around q = (0, 0) and q ~ (71", 0). 
We propose that the spin fluctuations for the d xy orbital 
around q — (0,0) and q ~ (tt,0) play an important role 
in the electronic states around x — 0.5 other than the 
criticality approaching the usual Mott transition where 
all electrons arc localized. 
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